GRB 011121: A collimated outflow into wind-blown surroundings 
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ABSTRACT 

We report optical and near-infrared follow-up observations of GRB 011121 collected predom- 
inantly at ESO telescopes in Chile. We discover a break in the afterglow light curve after 1.3 
days, which implies an initial jet opening angle of about 9°. The jet origin of this break is sup- 
ported by the fact that the spectral energy distribution is achromatic during the first four days. 
During later phases, GRB 011121 shows significant excess emission above the flux predicted by 
a power law, which we interpret as additional light from an underlying supernova. In particular, 
the spectral energy distribution of the optical transient approximately 2 weeks after the burst 
is clearly not of power-law type, but can be presented by a black body with a temperature of 
~6000 K. The deduced parameters for the decay slope as well as the spectral index favor a wind 
scenario, i.e. an outflow into a circum-burst environment shaped by the stellar wind of a massive 
GRB progenitor. Due to its low redshift of 2=0. 36, GRB 011121 has been the best example 
for the GRB-supernova connection until GRB 030329, and provides compelling evidence for a 
circum-burster wind region expected to exist if the progenitor was a massive star. 



Subject headings: gamma rays: bursts — techniques; photometric — supernovae: general 



1. Introduction 

In the presently favoured scenario for classical 
gamma-ray bursts (GRBs), occurring at cosmo- 
logical distances (measured redshift range so far 
0.36 < z < 4.5; van Paradijs etal. 1997; An- 
dersen etal. 2000), the explosion of a very mas- 
sive star leads to a fireball and a short, beamed 
flash of gamma-rays (Woosley 1993; Fryer et al. 
1999; Meszaros 2002), resulting in three physi- 
cally distinct observational phenomena, namely 
the GRB itself, a long-lasting afterglow and the 
classical supernova (SN) light. Whereas the af- 
terglow emission is probably fed by the kinetic 
energy of the coUimated, relativistic outflow, the 
supernova light is caused by the decay of radioac- 
tive nuclei created and released during the stellar 
explosion. The maximum of the supernova light 
is expected at ~10-20 {l+z) days after the ex- 
plosion, though at present it is not clear whether 
or not the GRB and the supernova explosion are 
delayed (Vietri & Stella 1998). One of the ba- 
sic consequences of this hypernova scenario is that 
it predicts (and links) the occurrence of a GRB 
with a jet, the unavoidable strong wind from the 
massive progenitor star, and the supernova light 
(Heger etal. 2003). 

Earlier evidence on the GRB-SN connection 
was initially based on the coincidence of GRB 
980425 and SN1998bw (Galama etal. 1998), and 
subsequently some supernova light contribution 
was found in the late-time light curves of GRB 
980326 (Castro-Tirado & Gorosabel 1999, Bloom 
etal. 1999), GRB 970228 (Galama etal. 2000), 
GRB 970508 (Sokolov etal. 2002) and possibly 
GRB 980703 (Holland etal. 2001), GRB 990712 
(Bjornsson etal. 2001), GRB 991208 (Castro- 
Tirado et al. 2001) and GRB 000911 (Lazzati et al. 
2001). 

Observations of massive stars in our Galaxy 
have shown that they lose matter via strong stellar 
winds in the LBV and in the Wolf-Rayet evolution- 
ary phase. GRB progenitors are expected to be 
very massive stars and perhaps physically related 
to Ib/c supernovae (Heger etal. 2003). The pro- 
genitors of these stars should develop fast winds 
at the end of their lives and thus, wind features 
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^ Based on observations collected at the European 
Southern Observatory, La Silla and Paranal, Chile (ESO 
Programme 165.H-0464). 



Fig. 1. — Light curve of GRB 011121 as recorded 
by HEND on Mars Odyssey in the outer scintilla- 
tor at >30 keV with 0.25 sec time resolution (top) 
and in the inner scintillator at >60 keV with 1 sec 
time resolution (bottom). The time axis zero point 
(UT) is indicated at the top. The three colours 
represent the three successive time segments used 
for the spectral fitting (see Fig. 2). 



are expected to be seen in afterglow light curves. 
The evolution of classical supernova remnants into 
wind-blown bubbles has been a target of detailed 
investigations since years (e.g., Benetti et al. 1999; 
Landecker et al. 1999). What has been missing 
in GRB research until GRB 030329 (Hjorth et al. 
2003), however, was a strong observational link 
between an underlying supernova component in 
an afterglow light curve and evidence for a fire- 
ball expanding into a stellar wind profile. Signif- 
icant work into this direction has been done over 
the past years by many groups from the theoret- 
ical as well as from the observational side (e.g., 
Meszaros, Rees, & Wijers 1998; Chevalier & Li 
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2000) and GRBs with wind-profile interactions 
have been identified in some cases, including GRB 
980425/SN 1998bw (Chevaher & Li 2000). How- 
ever, if one excludes GRB 980425 from this list, 
prior to GRB 011121 the observational data base 
was rather poor. The burst 011121 changed this 
observational situation, though in a less spectacu- 
lar way than GRB 030329. Not only was the dis- 
tance of the burster relatively small but also the 
light curve steepened 1 day after the burst so that 
a bright supernova component became visible. 

GRB 011121 was detected by the GRBM/WFC 
onboard BeppoSAX on 2001 November 21, 18:47:21 
UT and localized to initially 5 arcmin (Piro 
2001a). Subsequent analysis refined this posi- 
tion to 2 arcmin accuracy (Piro 2001b), and the 
triangulation of the GRB arrival times as mea- 
sured by Ulysses, BeppoSAX (Hurley etal. 2001) 
and HEND (Hurley et al. 2002) further refined the 
coordinates. A follow-up X-ray observation with 
the BeppoSAX narrow-field instruments revealed 
a fading X-ray afterglow (Piro etal. 2001). 

Follow-up optical/NIR observations of GRB 
011121 were quickly started by several groups, 
leading to independent discoveries of the after- 
glow (e.g., Wyrzykowski etal. 2001; Greiner etal. 

2001) . Further observations revealed a rather 
small distance, z « 0.36 (Infante etal. 2001), 
and excess emission above the early power-law de- 
cay (Garnavich etal. 2001). The interpretation of 
this excess light as a possible supernova bump at- 
tracted much attention and resulted in the so far 
best sampling of GRB afterglow emission at late 
times. 

Here we report the results of observations ob- 
tained by the GRACE^ (GRB Afterglow Collabo- 
ration at ESO) consortium. 

2. Observations 

2.1. X-rays/Gamma-rays 

The light curves of GRB 011121 measured with 
the two detectors of HEND are shown in Fig. 1. 
Integrated into the Third Interplanetary Network, 
the "Mars Odyssey" satellite has two instruments 
with GRB detection capabilities: the Gamma-Ray 
Spectrometer and the High Energy Neutron De- 
tector (HEND). HEND combines four detectors 

^http: //zon .wins .uva.nl/~grb/grace/ 



to measure the spectra of neutrons and gamma- 
rays (Mitrofanov etal. 2003). In particular, the 
inner/outer Csl-scintillators can measure gamma- 
ray photons in the 30-1000 keV range, at com- 
mandable integration times, and in case of a GRB- 
trigger (or solar flare), time histories at 1.0/0.25 
sec resolution are recorded. 

Fig. 1 shows that the gamma-ray emission is 
characterized by a main peak with ^10 sec du- 
ration, followed by a ^20 sec tail. Thus, GRB 
011121 clearly belongs to the long-duration sub- 
class of GRBs. 

1.1)00 p ' ' ' ' ' ' ' ' — 
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Fig. 2. — Energy spectra of GRB 011121 above 
200 keV, as recorded by HEND on Mars Odyssey 
in the outer scintillator for three successive time 
segments, shown as green, blue and red in the light 
curve (see Fig. 1.) 

According to the HEND energy spectra (Fig. 
2) there is no obvious spectral evolution along the 
light curve of GRB 011121. All three 20 sec seg- 
ments of the time profile have spectra with the 
same power law slope with a photon index of - 
2.35±0.25 (after backwards folding of the model 
with the instrument response). There is no evi- 
dence in the HEND data that the early afterglow 
emission started already at the time of the last 
segment of the burst. 

Based on the power law spectral slope of 2.35, 
the fluence is 2 ± 0.4 x 10~^ erg cm~^ in the 250- 
700 keV range. This corresponds to an (isotropic) 
energy release of = 2.7x10^^ erg at the given 
redshift (luminosity distance of 2.07 Gpc) and the 
cosmological parameters JIa = 0.7, — 0.3, and 
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Ho = 65 km s'^ Mpc'^ 





Fig. 3.— J band finding chart of GRB 011121 

obtained on 2001 Nov. 24 witli ISAAC/VLT. 
Marked arc the afterglow (AG), neighbouring ob- 
jects (1, 2), the host galaxy (3) and some local 
standards (J1-J4, B3, B8, B9). 

2.2. Optical and NIR Imaging 

About 9 hours after the GRB, Kg- and R- 
band images were obtained at La Silla (ESO) 
with SOFI (equipped with a 1024 x 1024 HgCdTe 
HAWAII array with 18.5 fim pixels and a plate 
scale of O'.'29/pixel) at the 3.58m NTT, and with 
EF0SC2 (equipped with a 2048 x 2048 thinned 
Loral/Lesser CCD with 15 /um pixels and a plate 
scale of O'.'157/pixel) at the 3.6m telescope, re- 
spectively. After our independent discovery of the 
afterglow, second epoch observations plus some 
exposures with additional filters were obtained 
(^ 4 hrs later) at the end of that night. Multi- 
colour photometry continued at the NTT and 
3.6 m telescope during the next three days and 
at the VLT thereafter. Due to the rapidly de- 
creasing afterglow brightness we switched to the 
instruments at the 8.2 m Very Large Telescopes 
(VLT) starting on day 3 after the GRB. For our 
near-infrared (NIR) observations, the short wave- 
length (0.9-2.5 /im) arm of the infrared spec- 
trometer ISAAC (equipped with a 1024x1024 



pixel Rockwell HgCdTc array and a plate scale 
of O'.'147/pixel) on the VLT- telescope Antu was 
used, and for the optical observations the Focal 
Reducer and Spectrograph FORS (equipped with 
a SITE CCD with 24 /zm pixels and a plate scale 
of O'.'2/pixel) at VLT/Yepun. 

Additional imaging was performed in the R 
and I bands at the ESO 1.54 m Danish tele- 
scope on La Silla (Chile), in the H band at 
the Anglo-Australian Telescope (AAT) and in 
the J and H bands at the 4-m Blanco tele- 
scope of the Cerro Tololo Interamerican Obser- 
vatory (CTIO). The 1.54m Danish telescope was 
equipped with a 2048x4096 EEV/MAT frame- 
transfer CCD (illuminated area is 2048x2048 
pixel) with 15 //m pixels, which provides a plate 
scale of O'.'39/pixel. At the AAT the IRIS2 in- 
strument with a 1024x1024 Rockwell HAWAII-1 
HgCdTe detector was used, having a pixel scale of 
O'.'446/pixel. At CTIO, the Ohio State InfraRed 
Imager/Spectrometer (OSIRIS) with a 1024x1024 
HAWAII-1 HgCdTc detector was used, having a 
plate scale of 0'.' 161 /pixel. 

2.3. Spectroscopy 

Spectroscopic observations were performed dur- 
ing the first two days in several wavelength bands 
to derive the redshift and during the maximum of 
the late-time bump to search for signatures of the 
potential supernova and to derive host galaxy pa- 
rameters. In the JK bands the short-wavelength 
arm of ISAAC was used at low resolution (with a 
final dispersion of 3.6 A/pixel in the J band, and 
7.2 A/pixel in the K band). The I'.'O slit was used. 
The JK band spectra were both taken at 1'.'2 see- 
ing, thus giving a FWHM resolution of 29 A (J 
band) and 59 A {K band). 

Optical spectra were taken with F0RS2 at the 
Yepun telescope at three different occasions, each 
with a different grism (see Tab. 1). The 1501 
grism has a mean 5.5A/pixel scale, which at the 
0'.'5 seeing and the use of a 0'.'7 slit led to a reso- 
lution of 14 A (FWHM). Grism 3001 has a mean 
2.59 A/pixel scale, leading to a resolution of 13 
A (FWHM) at 0'.'7-T.'2 seeing and use of the T.'O 
slit. Finally, grism 600B with a mean 1.2 A/pixel 
scale at 1'.'2-1'.'5 seeing and the I'.'O slit led to a 
FWHM resolution of 6 A. The pixel scale changes 
by less than 5% from the red to the blue end of 
each grism. The peak spectral response for the 
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Tabic 1: Log of the observations 



D.itr ;("Ti T.'l.' ■( ,)],(■ lii;-truujriit Tiltrr Cri-iu'"' E_viM.-ur- ;-.'<-i :-k'<-iii- D ri - lit m v- - i,nj;,g)(^) 



2001 


Nov 


22 


03:03- 


-03:25 


ESO NTT/SOFI 




10 X 63 


0''7 


15. 19±0.07 


2001 


Nov 


22 


04:03- 


-04: 14 


ESO 3.6m/EFOSC2 




3 X 180 


l''5 


18.86±0.02 


2001 


Nov 


22 


07:30- 


-08:01 


ESO 3.6 m/EFOSC2 


Rc 


3 X 600 


i"o 


19.51±0.01 


2001 


Nov 22 


07:47- 


-08:22 


ESO NTT/SOFI 


J 


25 X 70 


oi'8 


17.47±0.04 


2001 


Nov 22 


08:23- 


-08:47 


ESO NTT/SOFI 


H 


17X 70 


0"7 


16.72±0.05 


2001 


Nov 22 


08:47- 


-09:06 


ESO NTT/SOFI 


Ks 


13 X 70 


0"9 


16.09±0.07 


2001 


Nov 22 


14:00- 


-14:14 


AAO AAT/IRIS2 


H 


10X60 


2^'7 


17.30±0.25 


2001 


Nov 23 


00:04- 


-00:20 


ESO NTT/SOFI 




11 X 70 


l"2 


17.71ib0.10 


2001 


Nov 23 


00:20- 


-00:39 


ESO NTT/SOFI 


J 


13 X 70 


i^'e 


19.38±0.10 


2001 


Nov 23 


04:30- 


-04:40 


Hia\J 3.0 m/iiji' L>bO^ 


i?c 


3X180 


(c) 


21.23i0.05 


2001 


Nov 23 


04:41- 


-04:57 


ESO 3.6m/EFOSC2 


B 


3X300 


(c) 


22.40±0.50 


2001 


Nov 23 


04:58- 


-05:08 


ESO 3.6m/EFOSC2 


V 


3X180 


(c) 


21.95±0.15 


2001 


Nov 


23 


08:29- 


-08:45 


ESO NTT/SOFI 


J 


11X70 


0^'9 


19.73±0.10 


2001 


Nov 


23 


08:45- 


-09:01 


ESO NTT/SOFI 




11X70 


0"9 


18.29±0.10 


2001 


Nov 


23 


08:29- 


-08:39 


ESO 3.6m/EFOSC2 


Rc 


3X 180 


l"3 


21.75±n.08 


2001 


Nov 


23 


08:39- 


-08:44 


ESO 3.6m/EFOSC2 


B 


300 


l"6 


>21.70 


2001 


Nov 


23 


06:45- 


-08:34 


ESO Yepun/FORS2 


600B 


2x2400 + 1500 


l{'5,l"0,l"2 


— 


2001 


Nov 


23 


07:22- 


-08:34 


ESO Antu/ISAAC 


1.1-1.4 ^im 


3600 


l{'2 


— 


2001 


Nov 


24 


06:24- 


-06:48 


ESO Antu/ISAAC 


Ks 


15X60 


0^'7 


19.30±0.08 


2001 


Nov 24 06:53- 


-07:16 


ESO Antu/ISAAC 


H 


15X60 


0^'7 


19.81±0.05 


2001 


Nov 24 07:18- 


-07:36 


ESO Antu/ISAAC 


Js 


10X90 


o'.'e 


20.93±0.10 


2001 


Nov 24 07:53- 


-08:59 


ESO Antu/ISAAC 


1.95-2.55 /Ltm 


3600 


1. 2 




2001 


Nov 24 07:43- 


-08:01 


ESO Yepun/FORS2 


B 


3X300 


o'.'s 


23.70±0.30 


2001 


Nov 


24 


08:02- 


-08:14 


ESO Yepun/FORS2 


V 


3X180 


o'.'s 


23.37±0.05 


2001 


Nov 


24 


08:15- 


-08:27 


ESO Yepun/FORS2 


^special 


3x 180 


a"8 


22.93±0.08 


2001 


Nov 


25 


06:43- 


-06:59 


ESO Antu/ISAAC 


-/s 


10x90 


o'.'o 


21.64±0.30 


2001 


Nov 


25 


07:01- 


-07:24 


ESO Antu/ISAAC 


H 


15X60 


o'.'e 


20.73±0.10 


2001 


Nov 25 


07:27- 


-07:52 


ESO Antu/ISAAC 


Ks 


15X60 


o'.'s 


20.02±0.10 


2001 


Nov 


25 


06:34- 


-08:10 


ESO Yepuii/FORS2 


1501 


3 X 1800 


O'.'b 




2001 


Nov 


25 


08:10- 


-08:22 


ESO Ycpun/FORS2 


^special 


3X 180 


o'.'s 


23.68±0.15 


2001 


Nov 


25 


08:23- 


-08:34 


ESO Yepun/FORS2 


V 


3x 180 


(d) 


>24.40 


2001 


Nov 


25 


07:41- 


-08:13 


ESO 1.54 m Danish 


Rc 


1880 


l"l 


23.60±0.40 


2001 


Nov 


30 


08:05- 


-08:51 


CTIO OSIRIS 


H 


19x60 


o'.'i 


>19.40 


2001 


Dec 


01 


08:07- 


-08:39 


CTIO OSIRIS 


J 


8x120 


o'.'o 


>20.10 


2001 


Dec 


03 


07:50- 


-08:25 


CTIO OSIRIS 


J 


lOX 120 


l'.'2 


>20.10 


2001 


Dec 05 


06:28- 


-07:01 


ESO Melipal/FORSl 


V 


3X600 


0'.'9 


24.30±0.02 


2001 


Dec 05 


07:02- 


-07:33 


ESO Melipal/FORSl 


Rc 


3X600 


0'.'9 


23.20±0.0S 


2001 


Dec 05 


07:34- 


-08:07 


ESO Melipal/FORSl 


Tc 


3X600 


0'.'9 


22.35±0.25 


2001 


Dec 05 


06:53- 


-08:09 


ESO Antu/ISAAC 


Js 


40X90 


o'.'e 


22.41±0.15 


2001 


Dec 


07 


06:45- 


-08:02 


ESO Antu/ISAAC 


Js 


40X90 


0"6 


22.69±0.10 


2001 


Dec 


08 


06:20- 


-06:30 


ESO Melipal/FORSl 


V 


600 


o"s 


24.09±0.02 


2001 


Dec 


09 


06:49- 


-08:05 


ESO Antu/ISAAC 


Js 


40 X 90 


o'.'i 


22.52±0.15 


2001 


Dec 


09 


07:33- 


-07:43 


ESO A-Ielipal/FORSl 


V 


600 


i"n 


23.94±n.02 


2001 


Dec 


10 


06:47- 


-06:57 


ESO A-Ielipal/FORSl 


Rc 


600 


o'.'s 


23.65±n.08 


2001 


Dec 


11 


07:33- 


-08:06 


ESO Melipal/FORSl 


Rc 


3 X 600 


o'.'& 


23.45±0.08 


2001 


Dec 


11 


08:07- 


-08:29 


ESO Melipal/FORSl 


V 


2X600 


o'.'j 


24.18±0.02 


2001 


Dec 


12 


05:12- 


-07:30 


ESO Yepun/FORS2 


3001 


4X 1800 


o'.'i-\'.'i 




2001 


Dec 


13 


07:19- 


-08:46 


ESO Antu/ISAAC 


Js 


40X90 


O'.'l 


22.89±0.15 


2001 


Dec 


15 


05:10- 


-05:27 


ESO 1.54m Danish 


Ic 


4x900 


x'.'i 


>24.70 


2001 


Dec 


16 


04:51- 


-05:08 


ESO 1.54m Danish 


Ic 


2 x 900 


I'.'S 


>23.90 


2001 


Dec 


17 


04:55- 


-05:12 


ESO 1.54m Danish 


Ic 


4 x 900 


I'.'O 


>23.80 


2001 


Dec 


17 


07:20- 


-08:47 


ESO Antu/ISAAC 


Js 


40x90 


i"o 


22.79±0.30 


2001 


Dec 


18 


04:45- 


-05:01 


ESO 1.54m Danish 


Ic 


2x900 


l"l 


>23.80 


2002 


Feb 


09 


04:04- 


-05:24 


ESO Antu/ISAAC 


Js 


40x90 


0"5 


>24.85 



t*^) The -Rgpccial filter is about 10% broader than the standard Kc- Filters Js and are narrower with respect to the canonical J and K bands, respectively: 
Js has a width of 0.16 ^im (instead of 0.29 Mm), and has a width of 0.29 centered at 2.16 /Lim (instead of 0.35 jum centered at 2.20 /j-m). However, 

since both narrow-band filters have a higher transmission than the canonical filters, the net effect is that J — Js 0.05 (ftT — ftTs 0.02) when comparing 
the ISAAC Js vs. the SOFI J filter. The grisms 1501, 3001 and 600B are described in the text. 

Not corrected for Galactic foreground extinction, 
t"^) Image quality affected by guiding problems, 
t*^^ Variable seeing due to varying cirrus. 
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1501, 3001 and 600B grism is at 5000 A, 7500 A 
and 4000 A, respectively. 

Table 2; Local i)liotoinotrif D,V staudanls 



Star 


Coordinates (J2000.0) 


B 


V 


Bl 


ll'^34">18^9 


7G°01'37" 


20.11±0.01 


18.93±0.01 


B2 


ll''34'"19!7 


76°02'06" 


20.95±0.01 


19.95±0.02 


B3 


11^'34'"27!2 


76°01'52" 


22.93±0.03 


21.93±0.04 


B4 


ll''34'"27!2 


76°02'34" 


23.94±0.07 


22.49±0.05 


B5 


ll''34'"27!4 


76°02'22" 


21.58±0.01 


20.35±0.02 


B6 


lli'34'"29!9 


-76°02'36" 


zo.uyniu.uo 


99 97-i-n 0^ 


B7 


ll'>34'°33?3 


-76°01'33" 




1 q c;Q-i-n m 
icf.oymu.ui 


B8 


ll'^34'»36!2 


-76°01'07" 


zo.uumu.uo 


91 /T 1 _|_n OA 


B9 


ll''34'"38!9 


-76°01'38" 


90 '?7-l-0 ni 
zu.o ( inu.ui 


ly.uomu.ui 


BIO 


ll'>34'°41!9 


-76°01'04" 


9c: 9c:-|-0 1 


9Q 01 -i-n flQ 


Bll 
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Table 3: Local photometric Re- 


Ic standards 




Star 


Coordinates (J2000.0) 






Rl 


ll'^34"10!8 


-76°02'50" 


20.87±0.05 


20.06±0.02 


R2 


ll'^34'"ll!2 


-76°03'23" 


20.96±0.02 


20.35±0.02 


R3 


ll'>34'°15!4 


-76°00'37" 


19.97±0.04 


19.34±0.01 


R4 


ll'^34™18!0 


-76°02'32" 


23.43±0.10 


22.95±0.08 


R5 


ll'^34'»27!2 


76°02'34" 


21.75±0.04 


2L09±0.03 


R6 


ll'>34'»29!9 


-76°02'36" 


21.48±0.04 


20.80±0.03 


R7 


ll'^34'"33!3 


-76°02'30" 


21.82±0.06 


21.13±0.03 


R8 


ll'>34'°4ffi8 


-76°02'36" 


20.47±0.03 


19.98±0.01 


R9 


ll'>34'°41!9 


-76°01'35" 


22.36±0.05 


21.84±0.05 


RIO 


ll''34'"43!8 


-76°01'27" 


22.78±0.07 


21.49±0.04 


Rll 


ll'^34'"45M 


-76°01'32" 


19.36±0.08 


18.87±0.01 



R5=B4, R6=B6. 



3. Data reduction, analysis, and basic re- 
sults 

3.1. Photometry 

The optical and near-infrared images were re- 
duced in standard fashion using IRAF^ as well as 



^IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Associated Uni- 
versities for Researcii in Astronomy, Inc., under contract 
to the National Science Foundation. 



ESO's eclipse package (Dcvillard 2002). Photo- 
metric calibration of the GRB field was performed 
using SExtractor (Bertin & Arnouts 1996). Air- 
niass correction was done according to the coeffi- 
cients provided by ESO's Web pages (in units of 
mag airmass-i): ks = 0.240±0.007, ky = 0.112± 
0.005, kn = 0.091 ± 0.007, fc/ = 0.061 ± 0.006 for 
VLT FORSl and F0RS2, and fcj=0.06, fc//=0.06, 
A;if=0.07 for VLT ISAAC at Paranal, fcB=0.20, 
/cy=0.11, fcK=0.05, fc/=0.02 for La Silla. For the 
optical observations, the Landolt standard fields 
observed were SA 92-249, MarkA, Rubin 152, GD 
108 and SA 100. For the NIR calibration the 
UKIRT infrared standard star FS 12 was used. 

Local photometric standards (Tables 2, 3, 
4) were selected according to their detection in 
BVRI by SExtractor and a small scatter in their 
measured magnitudes on the images taken at dif- 
ferent observing epochs. Moreover, care was taken 
that the difference between the deduced Bessel-i? 
band magnitudes (FORSl) and the i?spcciai magni- 
tudes (F0RS2) of these stars was less than about 
0.1 mag. These local photometric standards were 
then used to derive the magnitude of the opti- 
cal transient (Tab. 1) after removal of the host 
galaxy. 

The accuracy of the photometry in the differ- 
ent filters has been cross-checked on the set of 
images taken on 24 November 2001 via a com- 
parison of the observed stellar brightnesses versus 
those of expected synthetic stars. For this we have 
first measured the magnitudes of our local stan- 
dard stars using YODA (Drory 2003), after having 
convolved the images to the common worst see- 
ing of 1'.'4. In a second step, we took synthetic 
stellar spectra from the Pickles spectral library 
and convolved those with the filter transmission 
curves and efficiencies of the corresponding instru- 
ment (F0RS2+ISAAC). Finally, we overplotted 
the measured standards over the synthetic stars 
in various color-color diagrams (Fig. 5). This 
shows that the photometric calibration is of co- 
herent quality over all the 7 filter bands used. 

3.2. The light curve 

GRB 011121 occurred at galactic coordi- 
nates (;ii,6ii) = 297?77,-12?43. The Schlegel, 
Finkbeiner, & Davis (1998) extinction maps pre- 
dict E{B — V) = 0.46 mag along this line of sight 
through the Galaxy. Assuming the usual ratio 
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Table 4: Local photometric standards for the NIR bands (Fig. 3) 



Star 


Coordinates (J2000.0) 




H 




Jl 


ll'^34™39!0 -76°01'17" 


15.21±0.05 


14.81±0.05 


14.91±0.05 


J2 


ll'^34"33!3 -76°01'33" 


17.28±0.05 


16.74±0.05 


16.82±0.05 


J3 


ll'^34"25!8 -76°01'37" 


18.56±0.05 


18.28±0.05 


18.27±0.05 


J4 


ll'^34"35!6 -76°01'53" 


17.15±0.05 


16.31±0.05 


16.31±0.05 




Fig. 4.— Sequence of J band images of GRB 011121 taken with SOFI/NTT on 22 November 2001 (left), 
and with VLT/ISAAC on 24 November 2001 (middle) and 9 February 2002 (right). The right-most image 
clearly shows the host galaxy (labeled #3 in Fig. 3), whereas the GRB afterglow has disappeared. The field 
is 17"xl7", North at the top and East to the left. 
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1 

3-V 



B-R 




Fig. 5. — Color-color diagrams of bright, non-saturated stars detected in the central area of the images of 
GRB 011121 (filled circles), compared with the expected color-color sequence as derived from the convolution 
of the stellar library of Pickles with the total filter efiiciency curves of the instruments used (open circles). 



of visual- to-selective extinction of 3.1, this gives 
Ay « 1.4 mag. The HI maps of Dickey & Lock- 
man (1990) give TVh = 1.2 x 10^^ cm'^ which 
translates into Ay « 0.9 mag using the extinction 
to absorption correlation of Predehl & Schmitt 
(1995). We used the higher of the above two val- 
ues, and then followed standard procedures (Rieke 
& Lebofsky 1985; Cardelli, Clayton, & Mathis 
1989; Reichart 2001) to calculate the extinction 
in the other photometric bands, resulting in: Ab 
= 1.87 mag, A^^ = 1.15 mag, Aj^ = 0.83 mag, 
Aj = 0.40 mag. Ah = 0.25 mag and Ar, = 0.17 
mag. In all the light curve plots this extinction 
has been corrected for. 

Accurate photometry of the afterglow requires 
a careful removal of contaminating light from the 



underlying GRB host galaxy which has an angular 
extent of two arcseconds (Fig. 4). This was ac- 
complished with a two-component intensity pro- 
file model, as described in more detail in §4.4. 
The model is derived from the J band image of 
9 February 2002, taken at a seeing of 0'.'45. We 
assume that the radial profile of the host galaxy is 
the same in all filter bands, as justified by a com- 
parison of the scale lengths of 86 face-on galaxies 
in the B,V,R,H filter bands by de Jong (1996), 
who finds that they are identical within the er- 
rors (the disk scale length ratio of R vs. H band 
is 1.07±0.08). The subtraction is done in count 
space on the images before the photometry, and 
for each image the host model is first convolved 
with the seeing of that exposure, and then itera- 
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tively adjusted in intensity until the area around 
the afterglow is best matched to the undisturbed 
background somewhat further away. As a consis- 
tency test wc computed the host magnitude for 
each image taking into account the seeing convo- 
lution and the final intensity normalization, and 
found that the host magnitude differs at most by 
0.07 mag (two instances) from the values given in 
§4.4. In addition, star 1 from Fig. 3 has been used 
to monitor the quality of the subtraction, which 
was always accurate to better than 0.05 mag. This 
ensures that the light curve shape is not contami- 
nated by this subtraction procedure. An example 
of the result is shown in Fig. 6. 

3.2.1. Light curve based on observations prior to 

t = 10 days 

Figure 7 shows the afterglow light curve in the 
i?c-band after removal of the flux from the under- 
lying host galaxy, and after correcting for galactic 
extinction. Theoretical fits were obtained with the 
formula of Beuermann et al. (1999) in the repre- 
sentation of Rhoads & Pruchter (2001): 

F,{t) = 2^/" F,{tb)[{t/tX' " + (i/*6)"' "]"'^" • 

(1) 

Here F,^ is the fiux density, tb the break time (in 
days), and n the parameter which describes the 
smoothness of the break. In all cases we fitted 
apparent magnitudes. Here and in the following 
we use the standard notation for the time and fre- 
quency dependence of the fiux density in the sim- 
ple fireball model: F^, oc t~°' y~l^ . 

The fit of the i?c-band data excluding the mea- 
surements at i > 10 days gives an early-time 
slope of ai = 1.62 ± 0.62, a late-time slope of 
a2 = 2.44 ± 0.38, and a break time tb = 1.2 ± 1.0 
days (x^/d.o.f. = 0.45; d.o.f. = degrees of free- 
dom) . The parameter n which measures the sharp- 
ness of the break, was fixed at n = 10, correspond- 
ing to a sharp break. However, values as low as n 
= 1 give similarly good fits. 

At early times (22 November 2001, 3:30 UT) 
we measure Rc — K = 2.7 ± 0.1 mag (after correc- 
tion for galactic extinction). Hence, the spectral 
slope of the afterglow is /3 = 0.80 ± 0.15. Two 
days later, on 24 November 2001, 7:00 UT, we 
find P = 0.62 ± 0.05 based on our BVRJHK 
data, which is consistent with the earlier result, 
and which is fully in agreement with the results 



obtained by Garnavich ct al. (2003) and Price ct al. 
(2002). These afterglow parameters are fully con- 
sistent with those expected from the simple ver- 
sions of the fireball as well as observations of 
previous gamma-ray burst afterglows (e.g., van 
Paradijs, Kouveliotou, & Wijers 2000). 

We can improve the fit by adding more data 
points as follows. First, using the deduced spec- 
tral index of the early afterglow we can estimate 
that at the time of the single AAO observations 
(Table 1), i.e. before the deduced break time, the 
i?c-band magnitude of the optical transient was 
20.12±0.25. Second, wc can include in our fit the 
i?-band data from Garnavich ct al. (2003) for 
t < 0.6 days. In doing so, we get ai = 1.62 ±0.39, 

= 2.44 ± 0.34, tb = 1.20 ± 0.75 days, xVd-o.f. 
= 1.11. 

3.2.2. Light curve based on all observations 

Including the "late-time" data (after t = 10 
days) to the fits adds substantially more freedom, 
since different model components can now com- 
pensate each other. As before, we only use Rc- 
band data for the analysis. In order to account for 
the excess light after t = 10 days, we assume that 
this "bump" is due to light from an underlying su- 
pernova. In fact, it has been officially designated 
as SN 2001ke (Garnavich etal. 2003). To model 
this supernova component we employ the observed 
UBVRJc light curves of SN 1998bw (Galama et 
al. 1998) as a template. 

The effects of the redshift of GRB 011121 
(^=0.36) were taken into account assuming the 
cosmological parameters as given in section 2.1. 
For SN 1998bw we used a redshift of ^ = 0.0085 
(Tinney etal. 1998). The entire numerical proce- 
dure is in detail explained in Zeh & Klose (2003). 
It interpolates smoothly between the UBVRI 
hght curves of SN 1998bw (Galama et al. 1998)^, 
so that for any given time a set of fiux densities is 
calculated. The predicted time- dependent super- 
nova light curve of a redshifted GRB-supernova 
based on the SN 1998bw template is then calcu- 
lated according to the procedure outlined by Dado 
etal. (2002a). The apparent magnitude of a red- 
shifted SN 1998bw in a given photometric band 
is finally obtained by integrating over the flux 
densities (transformed into per unit wavelength), 

*see also http://zon.wins.uva.nl/ titus/grb980425.html 
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Fig. 6. — Sequence of R band sub-images of 24 Nov. 2001 showing the original (left), the model of the host 
galaxy (center) and the host-subtracted image (right). All sub- images are shown at the same scale and with 
the same cut values. Labels are the same as in Fig. 3. Photometry of the lower-left star in the original and 
the host-subtracted image has been done for all epochs to ensure that no oversubtraction occurred. The 
residual scatter in this comparison star was in the range of 0.02-0.05 mag. 



multiplied by the corresponding filter response 
function, and by applying the usual normalization 
factors. The entire procedure has been succesfuUy 
applied for GRB 030329/SN 2003dh (Hjorth etal. 
2003; Zeh etal. 2003) and has been successfully 
tested against the results obtained by others for 
other bursts (e.g., Masetti etal. 2003; Dado etal. 
2002a). The results are shown in Figs. 8 and 9. 

The resulting parameters are: ai — 1.63 ± 
0.61, a2 = 2.73 ± 0.45, h = 1.26 ± 0.94 days and 
k = 0.85 ± 0.11 with x /d.o.f. = 1.11. Here, k is 
the luminosity ratio in the Rc band between the 
GRB 011121 supernova and SN 1998bw at max- 
imum (at z = 0.36). As before, the parameter 
n was held constant at 10 because otherwise the 
sparse data would not allow convergence of the 
numerical fitting algorithm. Our conclusions are 
not sensitive to the exact value of this parame- 
ter. On the other hand, we emphasize that the 
value we get for the parameter k assumes a 'per- 
fect' SN 1998bw light curve at the redshift of the 
burster. Any potential relation between luminos- 
ity and light curve shape of a supernova in a cer- 
tain photometric band is neglected here. In partic- 
ular, k does not say anything about the bolometric 
luminosity of the GRB supernova compared to SN 



1998bw. This parameter is also sensitive to poten- 
tial additional parameters which one can introduce 
in order to improve the fit of the SN light curve 
(see § 4.2.2). 

For SN 1998bw we assumed zero extinction 
along the line of sight through our Galaxy. Ac- 
cordingly, our results obtained for the luminosity 
of the supernova accompanying GRB 011121 in 
units of the luminosity of SN 1998bw, the parame- 
ter fc, scale ask^kx exp{-AviSN 1998bw)/1.086), 
if such an extinction is taken into account. Note 
also that we assumed a Galactic extinction along 
the line of sight of Ay(Gal)= 1.4 mag. If smaller 
values are preferred (Bloom etal. 2002), k has to 
be corrected/reduced again. 

Figures 4 (middle panel), 8 and 9 show that 
the SN light already affected the afterglow light a 
few days after the GRB. This implies a substantial 
impact on the measurement of the parameter a2 ■ 
Unfortunately, due to instrumental constraints, we 
have a data gap between days 4 and 10 after the 
burst trigger so that the la error of the deduced 
a2 rises substantially when we include the SN light 
in the fit. Note that most of the data published 
previously are not host subtracted, and therefore 
cannot be simply added to our data set (Brown 
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Fig. 7. — i?c-band light curve of the afterglow of 
GRB 011121 and a fit to the data prior to i = 10 
days. Compared to Tab. la conservative system- 
atic 1(7 error of 0.15 mag has been added to all 
data. The break in the power-law decay occurs at 
tt = 1-2 days. At late times the afterglow is signif- 
icantly brighter than predicted by the power-law. 
Note the particularly bright host galaxy. 

etal. 2001; Phillips etal. 2001; Price etal. 2001; 
Stanek & Wyrzykowski 2001). Despite this added 
tincertainty, our analysis clearly indicates the pres- 
ence of a break about one day after the burst. The 
exact time of the break is sensitive to the details 
of the fit (compare the results we obtain by fit- 
ting data prior to t = 10 days vs. fitting all the 
data), but the existence of the break at t ~ 1 day 
is clearly established. Previous studies of the af- 
terglow hght curve of GRB 011121 did not find 
this break (e.g., Price et al. 2002; Garnavich et 
al. 2003), which can be attributed to sparse sam- 
pling of the afterglow at this particular time. It 
remains difficult to explain why the radio data of 
the afterglow do not agree with a jet model (Price 
etal. 2002). 

3.3. The spectrum of the afterglow 

The VLT/F0RS2 spectrum of the afterglow 
(including light from the host galaxy and the un- 
derlying supernova) contains several strong emis- 
sion lines (Figs. 10, Tab. 5), but no absorption 
lines. The redshift determined from these lines is 
z = 0.362±0.001 (see the figures for hue identi- 
fications) , consistent with results of Infante et al. 



Fig. 8. — i?c-band hght curve based on eq. (1) 
and an additional supernova component, assuming 
dt = days (eq. 2). Note that all photometric 
data are corrected for Galactic extinction (§3.2). 
Compared to Tab. 1 a conservative systematic la 
error of 0.15 mag has been added to all data. 



(2001). This redshift corresponds to a luminosity 
distance of 2.07 Gpc (assuming the cosmological 
parameters of section 2.1), and a distance modulus 
of 41.59 mag. 

Using the foreground Ay (§3.2), we have de- 
termined the extinction values in the lines A^™°, 
and corrected the measured line fluxes according 
to i^corr = -^obs X exp(r), whcrc T ~ 1/1.086 x 
^iino |--|^g continuum extinction according to the 
galactic Ay. 

4. Results and Discussion 

4.1. The afterglow model and the a — /3 re- 
lations 

Since our finding of a jetted explosion is crucial 
for an understanding of the observational data we 
have carefully checked this result. When we fit our 
i?c-band data (Table 1) with a single power-law 
plus SN component we get a = 1.98 ± 0.11, fc = 
0.81 ± 0.11 with xVd-o.f. = 1.49. Not only is 
the fit worse in comparison to the fit based on 
the Beuermann equation, the value we now get 
for a is notably larger (3(t deviation) than the one 
deduced for the early-time slope of the afterglow 
light curve (Price et al. 2002; Garnavich et al. 
2003) . This again points to a change in the fading 
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Table 5: Measured line fluxes (corrected for galactic extinction) and luminosities. 



Line 


Aobs 


z 


Flux 


Luminosity 




(A) 




(10-1^ erg/cmVs) 


(10^° erg/s) 


Oil] 3726 


5079 


0.363 


0.41±0.03 


8.6 


H/3 


6625 


0.362 


0.22±0.04 / 0.14±0.04(i) 


3.2/2.0 


OIII] 4963 


6760 


0.362 


0.11±0.06 / 0.08±0.05(i) 


1.6/1.1 


OIII] 5007 


6821 


0.362 


0.35±0.05 / 0.16±0.04(i) 


4.9/2.2 


Ha 


8945 


0.363 


0.96±0.08 / 0.81±0.09(i) 


9.2/7.7 


Paa 


25550 


0.362 


2.58±0.12 


14.8 



The first number is measured from the 1501 spectrum on 25 November 2001, while the second number is 
measured from the 3001 spectrum on 12 December 2001. 




t (days after the burst) 



Fig. 9. — The same as Fig. 8 but now showing the 
best fit in all photometric bands. When fitting the 
V, J, H, Ks-h&nd light curves the afterglow param- 
eters deduced from the i?c-band fit (ai, a2, [n], tb) 
were used as input, i.e. the functional form of 
the afterglow light curve was fixed. The super- 
nova light curve (extrapolated towards the H and 
K bands), however, was calculated for the chosen 
photometric band. 

behavior of the afterglow between about 0.4 and 2 
days after the burst. Moreover, our deduced ai for 
the early-time slope of the afterglow light curve is 
fully consistent with the result obtained by Price 
et al. (2002) and Garnavich et al. (2003) based on 
their optical/NIR data. In fact, when we fit their 
data we confirm that they do not show evidence 
for a break in the light curve. We attribute this to 
the lack of data around t — I day. This makes us 
confident that the break in the light curve which 




4000 5000 
Wavelength (A) 



Fig. 10.— Spectrum of GRB 011121 in the blue 
band, taken on 23 Nov. 2001 with VLT/F0RS2 
equipped with the 600B grism. While only the 
[O II] line is detected beyond doubt, we also indi- 
cate the positions of other lines which could have 
been expected. 

we deduce based on our observations is a real ef- 
fect. 

In Tab. 6 we list the predicted spectral slope 
(3 corresponding to the measured a's for various 
afterglow scenarios. Basically, we have to decide 
here between an ISM and a wind model for a jet- 
ted explosion and to constrain the position of the 
cooling break frequency in the spectral energy dis- 
tribution (SED) of the afterglow light at the time 
of the observation. From the comparisons of the 
predicted /3, based on the theoretical a — f3 rela- 
tions, with the observed one on 22 November 2001, 
3:30 UT, and on 24 November 2001, 7:00 UT, we 
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conclude (see also Fig. 15): 

1. The model predictions agree best with the 

observations for the wind scenario. In this 
model, the presence of a late "bump" is the 
natural consequence of the supernova follow- 
ing the stage of rapid mass loss in a massive 
progenitor star. 

2. The data favor the interpretation that dur- 
ing our observations the cooling break fre- 
quency I'c, which separates the contribution 
of fast-cooling electrons from slow-cooling 
electrons in the afterglow light (see Sari et al. 
1998 for their original definition), was above 
the optical/NIR bands {uc ^ 10^^ Hz). This 
is in agreement with the conclusions drawn 
by Price et al. (2002) based on their observa- 
tional data; we refer the reader to this pub- 
lication for a further discussion of this point. 

3. The power-law index p of the electron en- 
ergy distribution, N{'yf.) d'ye oc 'y~P d'ye was 
close to the expectations from particle accel- 
eration in relativistic shocks (e.g.. Sari etal. 
1999). 

Within our measurement errors, there is no 
need for an additional extinction by dust in the 
GRB host galaxy. A wavelength-dependent ex- 
tinction by cosmic dust in the GRB host galaxy 
would tend to increase (3, i.e., redden the after- 
glow light. For example, this was observed in 
the optical/near-infrared afterglow of GRB 000418 
(Klose et al. 2000). The low value of P we and oth- 
ers (Garnavich etal. 2003: Price etal. 2002) find 
for the afterglow of GRB 011121 however, after 
correction for the influence of reddening by Galac- 
tic dust, compared with various model predictions 
based on the observed light curve shape (Tab. 6), 
gives us no strong hint for an additional dust 
component acting along the line of sight. Grey 
dust in the GRB host galaxy could still be there 
(wavelength-independent scattering cross section 
in the considered photometric bands), but this 
cannot be deduced from our data. 

On the other hand, one should be aware of the 
fact that the results obtained are sensitive to the 
adopted Galactic extinction. Ay {Gal), along the 
line of sight. For example, if we had used Ay {Gal) 
= 0.9 mag, as it suggested by the H I maps (§ 3.2), 



then the extinction corrected (3 for the intrinsic 
optical afterglow would have been 1.1, based on 
our measured Rc — K color of the optical tran- 
sient on 22 November 2001, 3:30 UT. We favor the 
higher extinction value (1.4 mag), however. The 
early multi-color UBVRIJK observations of the 
GRB afterglow by Garnavich et al. (2002) give a 
best fit for the Galactic reddening along the line 
of sight of E{B -V) = 0.43 ± 0.07 mag. For 
a standard ratio of total-to-selective extinction of 

3.1, this is basically consistent with our chosen 
value for Ay {Gal). We note also that an adopted 
extinction of Ay (Gal) = 1.4 mag is still consis- 
tent with the result deduced by Price et al. (2001) 
within their claimed Icr error bar {Av{Gal) = 1.16 
± 0.25 mag). 

We note that the observed change in the decay 
slope of the light curve around the break time is 
in agreement with the predictions of the jet-wind 
model, provided that the steepening of the light 
curve is due to the sideways expansion of the jet 
(Rhoads 1999). 

It has been argued (Kumar & Panaitescu 2000) 
that for an external density profile as the jet 
break in the light curve is expected to be very 
gradual, taking at least two decades in time before 
most of the steepening sets in. In the case of GRB 
011121, the break in the light curve is certainly 
shallower, extending no more than one decade in 
time (even with fits using n—1 the "break" extends 
only from 0.5-2.5 days, still less than a decade). 
While this may be considered a problem, more ex- 
tensive considerations have shown a large diver- 
sity of light curve shapes, depending not only on 
the density profile but also on the evolution of the 
Lorentz factor, whether jets are uniform or non- 
uniform, and the difference of viewing and jet an- 
gle (e.g., Wei & Jin 2003). 

4.2. The supernova and the wind signa- 
ture 

4-2.1. Supernova features 

The light curve of GRB 011121 provides the 
clearest case to date for excess emission above the 
usual power law extrapolation. The excess in this 
case appears to become significant about one week 
after the burst (Fig. 9). This feature has been 
heralded as strong evidence for a supernova com- 
ponent (Bloom et al. 2002; Dado et al. 2002b; Gar- 
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Table 6: Predicted (3 values for various afterglow scenarios. Assuming a relativistic jetted explosion then 

for observations aX t < ibrcak (pre-break time) the isotropic model holds and q = ai, whereas for t > tbrcak 
(post-break time) the jet model applies and a = a2- We use here ai = 1.62 ± 0.39 and 0:2 = 2.44 ± 0.34 
(see § 3.2.1). The parameter s is the power-law index of the density profile of the circumburst medium, 
n(r) oc r~^. For an ISM model s = 0, for a wind model s = 2. Case 1 stands for > i^c, case 2 for < i^c- 
In the former case the electron power-law index is given hy p = 2/3, whereas in the latter case p = 2P + 1 
(e.g.. Sari et al. 1999). 
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2.44±0.34 



navich etal. 2003; Price etal. 2002). Garnavich et 
al. (2003) proposed the supernova label SN 2001ke 
for this afterglow, although no obvious supernova 
features were apparent in their Magellan/LDSS2 
spectrum obtained on Dec 7, 2001. While it is 
far from proven that the late emission observed in 
GRB 011121 is in fact due to a supernova, this in- 
terpretation is most natural. If we believe the SN- 
bump picture, the observations of Garnavich et al. 
(2003) provide a stern warning: SNlQQSbw (used 
in our analysis) may not be appropriate as a "tem- 
plate" for SN-light associated with gamma-ray 
bursts. Given that to date we have a rather lim- 
ited sample of supernovae associated with GRBs 
(see the recent review by Weiler et al. 2002), we are 
limited in what we can conclude from the bump 
in the afterglow of GRB 011121, except that it is 
obvious that the simple afterglow model is unable 
to explain features like these. Shock rejuvenation 
and inliomogeneities in the GRB environment or 
energy injection can explain the kind of bumpy af- 
terglow observed in the case of GRB 021004 (e.g., 
Lazzati etal. 2002, Holland etal. 2003), but does 
not account for the long-term (weeks), sustained 
bump observed in GRB 011121. Supernovae, on 
the other hand, provide a natural explanation for 
the energy and also the time scale involved in the 
late bump we are considering here. 



During the maximum of the SN bump (Fig. 13), 
the observed brightnesses of the afterglow are J = 
22.8±0.1, R = 23.4±0.1, and V = 22.8±0.1. This 
is about 3.6 mag, 3.3 mag, and 3.7 mag fainter 
than the total light of the host galaxy (see be- 
low). Although not all of the host was covered 
by the slit, the host galaxy still dominated the 
collected light - even for the spectrum taken on 
12 December 2001, around the maximum of the 
bump. Thus, the lack of any supernova signature 
in the spectrum is no argument against the su- 
pernova interpretation of the bump (Fig. 11), as 
pointed out by Garnavich etal. (2003). 

The extinction-corrected (§3.2) absolute magni- 
tudes for the STipernova with a A'-correction in the 
VRJ bands of 0.85 mag (Leibmidgut 1990), 0.55 
mag and zero, respectively, are Mj = -19.18, Mr 
= -19.79, and My = - 19.83. Comparison with the 
absolute magnitudes of SN 1998bw suggest that 
the supernova in GRB 011121 is slightly brighter 
in the V band, but fainter at longer wavelengths. 

4-2.2. The GRB- supernova connection 

While we use SN 1998bw as a template super- 
nova, we added one degree of freedom to the fitting 
procedure by allowing for a shift dt of the observed 
SN maximum (^011121) with respect to the pre- 
dicted one for a redshifted SN 1998bw (Cggsbw)' 
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Fig. 11.— Spectrum of GRB 011121 in the red 
band, taken with VLT/F0RS2 equipped with the 
1501 (25 Nov. 2001, top) and 3001 (12 December 
2001, bottom; shifted by 0.01 units downwards) 
grisms, respectively. Some hne identifications are 
indicated. Taking sht losses on 12 December 2001 
into consideration, both spectra represent basi- 
cally the host spectrum, since at the covered wave- 
length range the afterglow has faded below the 
brightness of the host. The dashed line shows the 
spectrum of SN 1998bw 24 days after the maxi- 
mum (Patat etal. 2001), redshifted to z=0.36 and 
diminished in brightness to 85% (see §3.2.2), thus 
clearly demonstrating that the supernova was too 
faint with respect to the host galaxy to be discov- 
ered spectroscopically. 

i.e. 

II j.max #max /r\\ 

UL — ^011121 ''1998bw ■ 

With this extra degree of freedom added, the best 
fit parameters become: ai = 1.67 ± 0.58, a2 = 
3.55 ± 1.35, tfc = 1.34 ±0.60 days, /c = 0.86 ± 0.14, 
and dt = -6.7 ± 5.2 days with xVd-o.f. = 0.74. 
As before, the parameter n was fixed at n = 10 
(Figs. 12, 13). The fact that is less-well deter- 
mined when we include the data points obtained 
after day 10 is not surprising. It is basically due to 
the lack of data between day 4 and 10 combined 
with the early and rapidly dominating SN compo- 
nent that makes the error of the deduced a-i, rela- 
tively large. This is also the reason why in Tab. 6, 
where we discuss the appropriate afterglow model, 
we used the fitting results obtained for those data 
when the SN component is still negligible. 



Fig. 12.-— The same as Fig. 8 but with dt (eq. 2) 
as a further free parameter. 

The fit improves after inclusion of published 
late-time HST data (Bloom etal. 2002). With 
these additional data we find for the i?-band light 
curve: ai = 1.66 ±0.44,^2 = 3.43 ±0.68, 4 = 
1.33 ± 0.48 days, fc = 0.96 ± 0.06, and dt = 
-5.1 ± 1.85 days with xVd-o.f. = 0.60 (Figs. 14). 
Finally, if we follow section 3.2.1 and include the 
transformed _ff-band data point and the early-time 
data from Garnavich et al (2003) we get basically 
the same result with slightly reduced error bars 
for 0L\ and t\,. 

It was already previously noted (Bloom et al. 
2002) that a negative time delay dt (Eq. 2) of order 
of a few days provides a better fit to the data than 
setting = days, i.e. GRB and SN start at the 
same time. Basically, two scenarios could explain 
such a delay. First, this phenomenon could be in- 
trinsic to the ejected supernova shell itself in the 
sense that light curve shapes of type Ibc/II super- 
novae are a function of the mass of the progenitor 
and other details of the explosion (as already re- 
marked by Bloom et al. 2002). The more exciting 
alternative is that this delay could point to a gen- 
uine time delay between the supernova explosion 
(formation of a neutron star) and the GRB (in- 
terpreted as the subsequent formation of a black 
hole; see Vietri & Stella 2000). It is tempting to 
use the negative delay found in our fitting proce- 
dure as an argument in support of the supranova 
model (Vietri & Stella 2000). However, it is worth 
to remember that the supranova model requires 
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t (di^s after the hurst) 

Fig. 13.— The same as Fig. 9 but for dt = -6.7 
days. 

a delay between the supernova and the GRB of 
many weeks to months (Vietri & SteUa 2000), not 
days. If this delay is shortened, then the origi- 
nal goal of the model, namely to explain the long- 
lived iron lines, is dismissed. Much better data on 
both the GRB afterglow and the SN light curve 
are required to establish the reality of such an off- 
set. GRB 011121 has not been sampled enough 
to allow rigorous statements about the temporal 
relationship between SN and GRB, and the use of 
SN1998bw as a template is also a rather unreliable 
assumption of our modeling, as pointed out above, 
and further discussed below. 

4-2.3. Alternative explanations for extra light 

We note that there are three observational de- 
tails that do not fit the generally used template 
SN 1998bw: 

1. Rapid intensity decay: The J-band flux on 
9 February 2002 is far below the predic- 
tion of SN 1998bw. While there exist su- 
pernovae with rapidly decaying light curves, 
this faster decline applies only for times t > 
80 days. These supernovae are thought to 
possibly produce less Ni, but more Ti. In the 
case of GRB 011121, the observed faster de- 
cay happened within less than 60 days after 
the maximum, and thus cannot be explained 
with a low Ni production. 

2. Atypical color evolution: SN decay slower 
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Fig. 14. — The same as Fig. 12 but with the inclu- 
sion of published HST i?-band data from Bloom 
et al. (2002). 

at longer wavelength, e.g., SN 1998bw or 
SN 2002ap. The hght curve of GRB 011121 
shows that the brightness decay in the J 
band is very rapid. Moreover, the HST data 
demonstrate that the decay in the R band 
is not only slower than that in the J band, 
but also the F-band decay is slower than the 
i?-band decay. Thus, in GRB 011121 the 
color dependence is just inverted relative to 
a supernova: the longer the wavelength, the 
faster the decay. 

3. Different spectral energy distribution: The 
excess light in GRB 011121 is substantially 
bluer than that of SN 1998bw, as already 
noted by Garnavich ctal. (2003). This is 
most easily recognized in the different k 
values for the different filters: while k = 
0.86±0.14 for the Rc band (see previous sec- 
tion; excluding the HST data), we find k = 
1.10±0.50intheFband and A: = 0.40±0.03 
in the J band. 

We investigated the alternative to the super- 
nova interpretation, namely a dust echo due to 
scattered light (Esin & Blandford 2000; Waxman 
& Draine 2000; Reichart 2001). In the absence of 
extinction within the host galaxy, as in the case of 
GRB 011121, the dust echo is expected to be bluer 
than prescattered light, just as observed. This 
would suggest that a thermal dust echo is clearly 
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ruled out since it would peak in the near-infrared 
band (Reichart 2001). However, in the dust echo 
scenario there seems to be a severe problem with 
the temporal evolution. Dust echoes due to scat- 
tered light do not peak as sharply as thermal dust 
echoes or supernovae (Reichart 2001). The time 
dependence comes from the angular dependence 
of the escape probability of the scattered photons. 
Using the values of the differential escape prob- 
ability as presented by Esin & Blandford (2000), 
we find that the inferred temporal decay should be 
much smaller than the observed decay rate after 
the bump maximum. Other alternatives arc neu- 
trons in the blast wave (Beloborodov etal. 2003), 
the scenario of a refreshed relativistic shock in the 
GRB outflow caused by the wind profile of the 
GRB progenitor (Ramirez-Ruiz etal. 2001), the 
shell collision model (Kumar & Piran 2000) or con- 
tinuous injection (Bjornsson ct al. 2002). However, 
these scenarios are not yet detailed enough to be 
tested against observational data. 

4.3. The break in the light curve 

One might be concerned that the break in the 
R band light curve happens at an intensity level 
which just corresponds to the R band brightness 
of the host. However, this is chance coincidence, 
and only applies to the R band. There are three 
arguments in favor of the break being unaffected 
by the host brightness: (1) A subtraction of a con- 
stant flux could potentially produce a jump in the 
light curve, but no break with increasing devia- 
tion at later times. (2) The R band brightness 
of the host is the brightness integrated over the 
area of the host. In contrast, the additional flux 
which the host would add to the PSF area of the 
afterglow is at least a factor of 4 smaller. Thus, if 
the break would have been "produced" by a wrong 
host subtraction, it would occur at a level nearly 2 
magnitudes fainter. (3) The brightness of the host 
is different in different filter bands {V — J = 1.3 
mag) . This color is different than that of the after- 
glow {V — J = 2.45 mag; both values not corrected 
for extinction). Yet, the break is found at the same 
time in all filter bands, irrespective of the relative 
brightness of the host to the afterglow. 

The break time ti, we deduce has an la error 
of 0.5 days so that one might tentatively conclude 
that any break occurred between about 1 and 2 
days after the burst. Whereas a break at very 



early times (t < 1 day) is excluded by the data 
obtained by Price et al. (2002) and Garnavich et 
al. (2003), a break at t = 2 — 3 days could be hid- 
den by the bright SN component. The following 
discussion relies on the deduced break time of tb « 
1.3 days. 

According to Livio & Waxman (2000), for a jet- 
ted explosion into a wind-blown surroundings with 
an n(r) gas density profile it is 




where M_5 is the mass loss rate of the star in 
units of 10~^ Mq yr~^, i>w,3 is the wind velocity 
in units of 10^ km s~^ and S53 is the isotropic 
equivalent energy of the fireball in units of 10^^ 
erg. We introduced a factor (1 + z) in order to 
correct for the redshift. Using ^=0.36, £^53=0.27 
and assuming M-^/vvj.3, = 1 we obtain Ojot ~ 9°. 
This corresponds to a beaming factor (6 ~ 2/0?,^.) 
of about 100. 

Usually it is assumed that £'53 can be approx- 
imated by the isotropic equivalent energy release 
in the gamma-ray band during the burst phase. 
According to Garnavich et al. (2003), for GRB 
011121 it was = 0.27 x 10^^ gj.g (see also §3.4). 
After correcting for the beaming factor this gives 
an energy release of about 3 x 10"''" erg, in agree- 
ment with the typical energy releases of GRBs 
(PraU et al. 2001). 

Given the many multi-filter observations during 
the first 4 days (Tab. 1), we have constructed 
broad-band spectra from the optical to the near- 
infrared region. These are shown in Fig. 15 with 
a separation of about 24 hrs between each other. 
The spectral energy distributions of the first four 
days are consistent with a (3 = 0.70 power-law. 
This suggests that the break in the light curve was 
achromatic, in agreement with the predictions of 
the jet model. We note in passing that we do not 
see /3 to evolve towards an asymptotic value of 
1.1 as has been argued by Dado etal. (2002b) and 
assumed for the cannonball interpretation of the 
GRB 011121 afterglow light curve. 
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Nov 22 8:00 UT 




Fig. 15. — Extinction-corrected, broad-band spec- 
tral energy distribution of GRB 011121 on four 
consecutive days, starting about 14 hrs after the 
GRB, and on 5 December 2001 near the maxi- 
mum of the supernova-bump. The best-fit power 
law slopes (3 are given at the right side except for 
December 5. The maximum deviation of the data 
to the power law is ~0.15 mag, justifying the sys- 
tematic error added in the light curve fitting pro- 
cedure. It can be seen that the spectral slope does 
not change over the first four days, supporting the 
jet nature of the break seen in the light curve. The 
December 5 (filled circles) spectrum is clearly not 
of power law type. The smooth solid curve is a 
blackbody with a temperature of 6300 K, while 
the thin sohd line is the spectrum of SN 1998bw, 
redshifted to z=0.36. Note that on Nov. 24, ap- 
proximately 1 day after the break, the value of /3 
is well-defined and within the error bars well con- 
sistent with the predictions of the jet-wind model 
(Table 6; case 2). 

4.4. The host galaxy 

4-.4-1- Morphology 

Images in R, V, and J show that the host 
galaxy starts affecting the brightness estimate of 
the afterglow already after November 24, 2001. 
We have therefore taken the J-band image from 
February 9 to model the shape and intensity dis- 
tribution of the host galaxy in order to subtract 
the host flux from the earlier images. The Febru- 
ary 9 exposure consists of two sub-exposures, each 
with 30 min exposure time, and taken at a mean 
seeing of 0'.'45 and 0'.'55, respectively. We used the 



0'.'45 part and modelled the host galaxy using a 
bulge component 

^buigc(-R) = ^buigc.o exp 
and a disk component 



/ i? \i/4 



(4) 



Idisk{R) = /disk,o exp ( - 1.68 



(5) 



For each of these two components, the four param- 
eters /bulge, J (/disk.o)i effective radius Re, position 
angle (PA) and eccentricity (e = -\/l — (b/a)^) 
are derived using the IRAF/SPP package GIM2D 
(Simard 1998). The resulting fit parameters are: 

• Bulge component: /bulge, o = 14.78 mag/arcsec^, 
Re = 5.0 pixel, e = 0.19, PA = -41° 

• Disk component: /disk. o = 22.18 mag/arcsec^. 
Re = 19.5 pixel, e = 0.19, PA = -46° . 

The eccentricity derived from the galaxy profile in 
the J band corresponds to an inclination of ~10 
degrees. 

The fit to the radial profile is given in Fig. 16. 
To subtract the host galaxy emission from the 
images taken in November/December 2001, this 
model of the host galaxy was convolved with the 
corresponding seeing of each image, and then sub- 
tracted. For filters other than J, the total bright- 
ness of the model was modified to achieve good 
subtraction in the corresponding band, i.e. the 
normalization was varied to minimize the residu- 
als. Aperture photometry on the subtracted im- 
ages was then applied to obtain the magnitudes of 
the afterglow (see Tab. 1 and Fig. 9). 

Using the model parameters given above, and 
integrating over the galaxy, we obtain the follow- 
ing apparent magnitudes of the host galaxy: J = 
19.1±0.1, R = 20.1±0.1, and V = 20.4±0.1. 

The decomposition of the surface brightness 
profile of the galaxy into a disk and a bulge com- 
ponent shows that the bulge dominates within the 
inner 0'.'5 only. The bulge contributes about 22% 
to the total light. This can be interpreted as a 
Sbc morphological type. The observed bright- 
nesses of the host galaxy in the V, R and J 
band correspond, after correction for foreground 
absorption and applying the fc-correction, to ab- 
solute magnitudes of -22.58±0.53, -22.24±0.24, 
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Fig. 16. — Surface brightness profile of the host 

galaxy of GRB 011121 in the J band, measured 
on the 2002 February 9 image (first part) at 0'.'45 
seeing with a plate scale of O'.'147/pixel (two data 
points per pixel). The abscissa is the distance 
from the center measured as \J a? + 6^, with a and 
h being the major and minor axes of the galaxy. 
The dashed line denotes the bulge component, the 
long-dashed line the disk component, and the full 
line the sum of both components. 

and -22.09±0.25, respectively. This is at the 
very bright end of the range populated by spirals 
(e.g., Binggeli et al. 1988); the host galaxy of GRB 
011121 is thus among the brightest 5% Sbc galax- 
ies. 

The gamma-ray burst happened at an offset of 
0'.'9 from the center of the host galaxy (as mea- 
sured on a VLT image taken at 0'.'6 seeing), corre- 
sponding to a physical distance of 9 kpc (see also 
Ryder et al. 2001). The above derived disk radius 
is i?e = 29 kpc. This offset location of the gamma- 
ray burst with its supernova is consistent with the 
general supernova picture, as well as with offsets 
found for earlier GRBs (Bloom etal. 2002). 

The star-formation rate 

Using the extinction corrected line fluxes of 
[O II] and Ha (Tab. 5) and the relations given by 
Kennicutt (1998), i.e. SFR (M/yr) = 1.4x10-*^ 
X L(OII) and SFR (M/yr) = 7.9x10-^2 ^ L(Ha), 
we deduce a star formation rate of 1.2 M0/yr and 
0.72/0.61 Mo/yr, respectively. In contrast, Sub- 
rahmanyan etal. (2001) derive a star-formation 
rate of 13-70 Mq /yr from their radio detection of 



the host galaxy at 0.05 mJy and assuming a spec- 
tral slope of -0.5 between 1.4 and 8.5 GHz. The 
rate derived here is at the low end of the range 
found previously in GRB host galaxies. However, 
also for these earlier cases different results are ob- 
tained for different methods/lines, in particular 
also a large difference between radio- vs. optically- 
based data (Vreeswijk etal. 2001). 

We have no explanation for the fact that the 
host galaxy appears to be among the brightest 5% 
Sbc galaxies, and yet the derived star-formation 
rate is rather low. Also, the lack of clear evidence 
for dust in the GRB host along the line of sight 
is somewhat surprising in the supernova picture, 
though there is also no strong intrinsic extinction 
in most previous GRB hosts. 

5. Concluding remarks 

While our a and P values are consistent within 
the errors with those reported by Garnavich et al. 
(2003) and Price etal. (2002), the change in our 
values due to the discovery of the break and the 
evolution during the four days after the GRB leads 
to a preference of the jet-wind model. The ob- 
served break and its interpretation is in contra- 
diction with the interpretation of the radio data 
by Price et al. (2002) suggesting a break at times 
later than 8 days after the GRB. We cannot offer 
a solution to this problem, as the break is observa- 
tionally evident from a dozen of data points spread 
over several filter bands. 

Our observations reveal a light curve break at 
early times and the appearance of extra light at 
late times. The former is believed to be evidence 
for a collimated outflow (a fireball jet), whereas 
the latter could be due to supernova light. Fi- 
nally, the a — /3-relations favor a wind model 
(Table 6). The afterglow of GRB 011121 thus 
provides another case in support of the current 
standard paradigm of long-duration gamma-ray 
bursts: the collapsar model (Woosley 1993, Hart- 
mann & Woosley 1995) in which the birth of a 
black hole inside a rapidly rotating massive star is 
announced via jet formation, breakout, and propa- 
gation into a stellar environment that was shaped 
by the strong winds of the pre-coUapse star. In 
fact, GRB 011121 is the first case, where all three 
signatures (jet break, wind density profile and su- 
pernova bump) have been found, while in earlier 
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cases only two of these three ingredients could be 
found (e.g., Jaunsen etal. 2001). However, even 
with GRB 011121 this line of reasoning remains 
circumstantial, and many more afterglows have to 
be observed, and better sampled than is presently 
the case, to seriously establish the GRB-coUapsar 
picture. 
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